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i. INTRODUCTION 


In homogeneous reactor theory the thermal neutron flux 
distribution in a given direction in tne reactor core 18 
represented as a smooth curve havinz a shape which is de- 
pendent upon the core geometry., In a bare rectangular 
parallelepinod reactor the flux distribution along each co- 
ordinate axis Is a cosine curve acmss the core, Ina 
heterogeneous reactor consisting of rods of fuel remolarly 
arransed in a moderating medium, the ceneral shape of the 
flux distribution across the core is similar to the cosine 
distribution of the homoreneous reactor. However, the abe 
sorption of neutrons is much hisher in the uranium fuel rods 
than it is in the moderating material, and hence there are 
local depressions in the neutron flux near the uranium fuel 
rods. 

The theory of tone natural uranium neterogeneous réeuctor 
has been broken down into microscopic theory and macroscopic 
theory. Macroscopic theory deals with the overall flux äis- 
tribution in the reactor and permite the determination of 
such parameters as critical sige and critical masa for a 
given reactor design. Microscopic theory deals with the 
local flux distribution in the unit cell of the reactor core, 


and it permits the calculation of the various lattice 
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constants, such as thermal utilization, resonance escape 
probability, lattice diffusion lengtn and material buckling. 
The subcritical assembly can be used to determine experi- 
mentally these lattice constants for a proposed reactor de- 
sign. “ince reactor theory is subject to many limitations 
and approximations, the subcritical assembly is a valuable 
tool which can be used to either check or supplement theo. 
retical calculations. 

The purpose of this thesis was to investigate the flux 
distribution in the unit cell of the Towa State College 
uranium craphite subcritical assembly. ‘Several ‘techniques 
for flux measurement using the foii activation method were 
also investizated, The flux distribution was measured in 
three different directions inside the unit cell both with 
and without coolant, and the experimental results were com- 
pared with the theoretical flux distribution in the unit 
cell, 
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II. REVIY OF TIN LITARATURE 


The use of the activation method for measuring thermal 
neutron flux was covered In detail by Feld (h), Cohen (3) 
described further use of the activation method and how it 
could be used to determine the miecroscopie flux distribution 
within a unit cell. He also pointed out the particular 
value of the foil activation method in determining flux disg- 
tribution near the boundary of two dissimilar mediums where 
diffusion theory cannot be applied with accuracy. iiunmel 
and Kamerasah (9) investigated the neutron flux depression 
in the neighborhood of a silver foil but no quantitative re- 
sults were given for indium foil. Clayton (2) and Richey 
(11) discussed in detail the foil placement in the unit cell 
of a uranium graphite lattice and the procedures and correc- 
tions used in counting. 

Murray (10) developed the flux distribution for a two- 
region fuel-moderator lattice systen based on diffusion 
theory for monoenerretic neutrons. He further presented a 
method of estimating the effect of extra absorption due to 
the presence of other components in the unit cell, such as 
Gladding, tubins, coolant and insulation. In this method it 
was assumed that all the other components act as poisons which 


can be tolerated, and nence they do not appreciably disturb 
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The thermal neutron flux in a fairly lirge subcritical 
gareroly in the central regi n away from ti. boundaries and 


extraneous neutron sources can ve represented by (6, p. 261) 
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where f is tha thermal neutron flux, WV is the Laplacian 


Ec ` 
operator and D, Is the material mekling of theo particular 
inttice systen, usually oxpressed in om "2, vita the usual 
boundary conditions that the flux is everywiere finite and 


non-negative and is zero at bhe extrapolated boundaries tne 


solution to tha above eountion is 
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ae Two „egion System 


A first avoproximetion for the thermal neutron flux in 
tne unit cell can De made oy use of one proun diffusion 
theory in a two-region fuel-moderator system (10). To 
simplify the mathenstics the square cell is replaced by a 
cylindrical cell of equal area. The dúiffusion equatíon for 


tne fuel is 
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neutron flux in noutrons/om^ sec, 2. is tne macroscopío 

aosorotíion cross section in mt, and 3 is the thermal 


neutron oroduction rate ver om^. With the assumptions that 
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Sy = Ü, Sy is constant, AZ does not vary along the cell axis, 


and that Z ls constant at any given coll radius, the solu- 
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In the above ecuastions X is the inverse diffusion length for 


tho given medium and {ορ 14; K and X. are modified vessel 
functions of the zero and first order. The physical con- 

stanta for the suberitical assemvly were previously deter- 
mined (8) and are listed in table 8 along with the various 


1 


dimensions of the unit cell,  .ith tnese constants and a 
table of Bessel functions (1) the flux in the fuel and in 


the moderator was determined fron ይባ8» li and le assuring 
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neutro | | 
y = 1 πηγε . The theoretical flux distribution 


em“ ne G 


normalised to toe flux at tie cell boundary is shown in 


Figure 31. 
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The theoretical Flux üirtribution in a multirerfion unit 
cell consistin;: of usl, cladding, water coolant, aluminum 
process zubo and Frapäalte moderator is seased on the thermal 


utilization equation Jevived by Aumsey and Velkoff (5, 12), 
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Figure l. The sudcritical ass 
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placed on each side of the source. The center tank con- 
sisted of three comoartmenta. The two end compartments, each 
about 26 in. long, were filled with water. The center com- 
partment was left dıy, and in it was placed a small table on 


which the sources were mounted. 
η €> 
& e 0uli' ces 


The assecvbly source consisted of five individual vlu- 
tonium-beryllium neutron sources, each emitting approximately 
4463 ሼ 10° neutrons per second. Lach source was contained 
in a "teinless steel and tantalum container which was one 
inch in dlametor aná 1 3/8 in. high. ‘nen placed on the 
small source table which was located undernesth the center of 
the asgemoly the tops of the source containers were about 
1/16 in. beneath the floor of the assembly. The five sources 
were arransed in a cruciform shape oriented on the x and y 


axes Of the coordinate system used ag shown in “igure 2. 
3. fuel elements 


1he ansev»bly was loaded wito fuel elements as shown in 
Figure l by filling every other hole giving an 8.l8-1In, 
2quare Lattice in the Lower region of thm assembly. The 
fuel assembly consisted of canned natural uranium slugs 
wrapped with 2S aluminum wire spacers and inse:ted in 615 


aluminum process tubos. The uranium fuel itself consisted 
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of rods 1 in. in diameter and ὃ in. long. Tho 23 aluminum 
earns had a 0.040-in, well thickness and the end capa were 
0.200 in. thick. Thus the overall dimensions of the fuel 
sluus were 8.40 in. long by 1.060 in. in dianeter. Seven 
slugs were inserted in each process tube, The aluminum 
process tuoes were 62 in. long, had an outside diameter of 
1.375 in. and a wall thickness of 0.035 in. The effective 
thickness of the coolant annulus between the slug and process 
tube waa 0.112 in. The aluminum wire spacer was 0.10% in, 
in diameter, and approximately ten feet of wire was used in 
each fuel assembly. The onds of tre process tutes were 
plugged with number seven rubber stoppers when making runs 


with coolant. 
h. Indium foil positions 


Slots for inserting indium foils were located as shown 
in Figure 2, anc were used in making vertical and herizontal 
flux surveys of the overall assembly. A grid system vas used 
in identifying blocks and/or foil positions in the assembly. 
Layera were mimbered from vottom to too from 1 to lh and the 
ten columns were designated A through J from left te right on 
the east face of the assembly. The folis normally used for 
pile surveys were 1.0 inch by 1.5 in. by 0.003 in. thiskx and 
veizhed approximately 0.5 mg each. These were mounted on 


aluninus backing and vere inserted in the cile by msans of en 
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aluminum strin ቿጅፓ።31 ከ01 ;5፲ጀ5 It ws thus nossiule to ootain 
guPvoys At x = -3 ino, 2 7 30 in. end ut any value of y 


setueon zero and -309 in. 


er 


lock D-6 was cut vertically at a point 20 in. ín fron 
the east face of the assembly to provide a test section at 
teh unit cell flux seasurements could oe wade. This 
particular position was selocted to keep harmonic eifects to 
a minimum. The blecirs avove bloey D-6 were supported by a 
lever arrangement so buat one third of vlock ved could easily 


m moved in and out of the pile, Urooves in. deep and 


ዳጋ 


0.015 in. wide wore cub into tne saved-off face of tne 


graphite block spaced 3/4 in. apart. The unit oell turetker 


δα, 


with the foil positions is shown in ®igire 3. un the T 

A Phdials bhese were seven foil nositionr within the fuel 
aasemoly numsered fron 1 to 7 as shown for tho r radial in 
elmre 3. Since there was no gir space between tho »roooss 
tube and the graphite on tie . radial tnere were only four 
foil positio nurbered .1 through ca, within the ruel 
aggenbóly on this radial. The foil positions in tne grapalte 
were Maveced con ecutively nrocsedäing out the resusctive 
racdiel as shown in Figure 3. “oil positions along tie aid 


L. Padials extended to the unit cel) oouoóonury mile those 
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Figure 3. Unit cell 





along the 4% radial extended to tre fuel assembly in the ad- 
jacent cell. The coordinates of each foil position are 
listed in Table 1. A process tube was cut et a point 20 in. 
in from the east face to nermit olacing foils inside the 


coolant annulus of the fuel assembly. Un runs made witn 





table 1. Unit cell foll positions and vosition correction 
factors 
Position x z nadial fx Pa 
divtance 
in. in ins 
0.0 -12.00 230.00 ዕ 1.060 1.000 
Pl “12.00 30.54 0.54 1,080 1.034 
ጅ2 -12,00 30.60 0.60 1.000 leOhB 
P3 -12.00 30.65 0565 1.060 1.052 
nb -12.00 30.69 0.69 1.000 14 052 
P5 -12.00 30.75 0.75 3.009 1.052 
P6 -12,00 30.91 0.81 1.000 1⁄4962 
FT 212,00 30.87 0.67 1.000 1.066 
፣፡8 -12,00 31.50 1.50 1,080 1.119 
P11 -12.00 33.75 3.75 1.000 1.300 
Ρ12 -12,00 34.50 11,50 1.000 1.360 
al -11.62 30, 38 0.54 0.985 1.026 
qe -11.58 30.12 0.50 0.985 1.030 
48 -11.54 30.46 0.65 0.985 1.030 
wel -11,51 30.59 0.59 9.982 1.034 
OoOo HR 30175 1.06 0.97% 1,087 
36 “1050 31.50 2.32 0.552 1.110 
«T & 9.05 32.25 3.18 04533 Lele 
a9 Be 7.7. 5.30 0.898 1.300 
¿10 - 7.50 234.50 6.36 0,882 1.360 




















Table 1. (vontimued) 
sosition x š hadial Pa “ሐ 
distance = 
in. in, on. 
Q11 - 6.75 35.25 7 52 0.870 1.440 
E “11.46 30 «Ὁ "Ln 0.380 1.000 
ne -11.40 30,0 0.60 08.7793 1.000 
a = 11.5 33. 30.0 0.65 0.975 1,600 
¿5 ~11.25 30.0 0.75 0.972 1.000 
nb =1 1 19 30.9 0.61 0.971 1.000 
= -11.13 30.0 0. 7? 0.968 L. 030 
ና ~ 9.75 30.0 2.25 0.931 1,000 
(14.1. - 0.25 30.0 3. 75 0.998 1.000 
nio - 7.50 30.0 11.50 0.082 1, 900 
R13 - 6.75 30.00 5.25 0.879 1.000 


coolant the process tube vas ¡sealed with waternool alec- 
trician'!s tape. 
Three 53268 of Indiun foil were used for flux measure- 


ments in the unit cell as follows; 


Table 2. Indium foils used in unit cell 

pol 3458 (in.) Average wt. (mg) 
3mall Ë X Š 0.096 
Medium E x 3/4 0,130 
Large EX 7/8 01152 








“ት 
Ciy 


The weight of each foil was determined to t.e nearort tenth 
of a milligram. The foils were mounted on scotch tape 
backing and were held in place in positions «around the fuel 
element by means of electrician's tape or adhesive tape. 
kadial positions 2, 5, 6 and 7 vere obtained by Sending the 
tape into an inverted "U" with the foil placed at the 
desired position. 


Ge Counting Equipment 


£ 


À Nuclear-Chbicapgo model 161 A scaler and model 23) mica 
end window counter were used to count irradiated indium foil 
activities. “he counter was placed inside & 2-1n. thiok 
lead shield which resulted in am averaie cackground count of 
20 counts per minute. An automatic timer which could be set 
for any desired counting time was used in conjunction with 
the scaler. The indium foil counting geometry vas held 


constant by means of trays on which the foil positions had 


been marked. 





VI. LAÄPEUIMENTAL PIotch Alk 


A. Determination of / 


in order to correct foil readings obtained at various 
pointa in the suberitical assembly, it was necessary to 
determine 7 , the inverse relaxation length for the thermal 
neutron flux in the assenbly. Vertical flux surveys were 
made at x = -3 in. and y = -10 in. from a = 18 in. to 2 = 
Sk in. Points for z less thon 10 in. and greater than 5h 
in. were not used due to the proximity of the source jn the 
first instance and the change in lattice size in the latter. 
The indium foils weichiny an average Of 0.5953 gm vere used 
for these surveys, and they were irradiated for a minimum of 
eight hours which gave an induced activity of 99.6 per cant 
of the saturation activity. Observed activities were cor- 
rected back to time of removal from the assembly, and this 
saturation activity was then divided by the partícular foil 
weight to give the normalised saturation activity, 4,,, in 
counts per minute per gram of indium, ‘Surveys were made with 
and without water in the coolant annuli. -vounting times were 
adjusted to keep the relative standard deviation of the 
observed counting rate less than 4 per cent. 


The normalized saturation activities were plotted on 





seml-logarithmie paper and a straight line was faired 
through the points. The slove of this line yielded a tria 
value τον Υ 33, which was now used to compute the hermonic 
and end correction terms, Ce and ጭ* These correction terns 
were then divided into the normalized saturation activities 
to give Âq]? the activities which would ve obtained if the 
l,l harmonic of the flux distribution were the only one 
present. “o further refine the value a it was neces- 
sary to uso an iterative process whereby new correction terms 
would 355 computed ma epplied to the orluinsi normalized 
saturation acetiviti>s to obtain new corrected values of 511 
Tho method of least squares was applied to obtein s new 
value for πε. ror the purposes of this investigation suf- 
icient accuracy in the value cf ጊጊ gas obtained by going 
through the iterative procedure only once. Harmonic effects 
beyond the third harmonic were Pound to be negligible and 
were ignored in caleulating the harmonic — terms. 
Sinilariy the end corraction tern, Ga» was found to have 


magilgible effect beyond the first narmonic, so that G 
-2 Yaa 6-3) 


o Was 


assumed to 56 simply lee where c was taken esual 
to 79 in., tne heisht of the assembly, 

The values of wa Yor the higher harmonics were calcu} 
lated fron the relation (2) 


L 


L 
77 de 
Tam les (£) (m + n -= 2) + J Zac. 33 





27 


Table 3. inverse relaxation lengta ard buckling 





























without coolant With ecolant 
ፖር (in^) 0,0705 0.0713 
h3* 75 (in™*) 0.1598 0.1600 
733 (587*) 0,215 0,235 
Be (in”“) 1.0 x 10% 0.6 x 1074 
μέ (em"^) 28.8 x 107° aa y" 








ope caen sif 





here πα is the length of the side of the svuareebased 
ns5enbly including the extrapolation distance. The value of 
A vas measured to be 62 in. The values of Yan for the 
various hermonics with and without coolant are listed in 
Table 3 together with the values for the buckling. The 
material buckline was evaluated from the ecauation for a 
Square=based assembly, 


re 
መ 
Lo 


2 
ደ TT Í ኀ፤ 
5: “(213 - Λι ቁጥ ጋ 


so 


The values of harmonic and ond correction terms, La and 
Cys The normalized saturation activity, > w: and the cor- 
rected activity due to the first mode only, Ay], are listed 
in Table L, The corrected activities, ጳገጊ» Are plotted in 


Figures h and 5. 
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Tablo h. Vertícal flux survey at x = =3 in., È -10 in. 
Fosition g “6 Ch Ce En Á o ALT 
(in.) (c/m) ( c/m) 
without coolant 
5ን 16 0.9997 1.0047 1.0847 2750 2540 
14. 24 0.9991 1.0491 1.0491 19h 1851 
፲5 30 0.9989 1.0265 1.0280 108 1051 
E6 35 0,9974 1,0168 1.013 She 831 
E7 112 0,9911 1.0097 1, 902 535 534 
18 148 0.9810 1.0057 0.987 339 313 
E9 sh 0, 2680 1.0037 0.4972 101 158 
with coolant 
53 18 0.9990 1.0869 1.0869 2860 2635 
th 2h 0.9996 1.0505 1.0505 1862 1772 
ES 30 0› 9959 1.0296 1.029 1153 1121 
EG 36 0› 7975 1.0171 1.013 750 110 
R7 12 0.9945 1.0109 1,991. 503 501 
1.8 LB 0.9870 1.0058 0.991, 306 310 
2 5 0.9700 1,0034 0,971. 193 196 


B. Correction Factors for 





Unit Cell Poil Positions 


Harmonica and end correction factors, fa anda P S? wore 
calculated for exch foil position in tho unit celi. It 
should be noted that the correction factor ís equal to the 


reeinrocal of tie correction term 


fo = ኗ j fc u Eq. 35 


Harmonica snd end correction factors for each foil position in 











Pigure h. 


Vertical flux survey without coolant 


The vertical »ils surve; was made ab x = =3 in., 
y = -10 in, with 1 in. by 1s in. foils. The 
vertical unit cell survey was made at x = -10 
in., Y = -10 in, with & in. by 3/4 in. foils 

&t spacing 2. Unit cell survey data was 
normalized to pile survey data. 
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Pigure 5. Vertical flux survey with water coolant 


ME 


The vertical pile survey was mado at x = -3 in., 
y = -10 in. with 1 in. by 1% in. foils. The 
vertical unit cell survey was made at x = -10 
ine, y = -10 in. with $ in. by 3/4 in. foils 

at spacing 2. Unit cell survey data was 
normalized to pile survey data. 
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the unit cell aro listed in Table 5. 
in order to compare the flux distrivoution in tno unit 

cell with the theoretical flux distribution it was necessary 
to convert she activities at the various foil positions to a 

common reference point. The point chosen was tue center of 
the unit cell sxanined which corresponded to the center of 
the uranium slug at x * -12 in., z = 30 in. It was therefore 
necessary to sake corrections to sil activities or the co- 
Sine distribution in the x direction and for the exponential 
decrease ín the 2 direction. These position corrections were 


called f, and f, respectively, and they were evaluated from 
el 


Tx lex 
cos & κ--12 (35 


fz = πχ o >= 
cog UR cos 
xx 


t. 86 HG. 37 


the equations 





UN 


The above factors were muiltiplied together to pive one over- 
all correction fnctor, i, for each position in the unit cell 


&s follows 


Values of the position correctian factors are listed in aule 


7 


1, end values of F are listed in Table 5, 





Table 5 5 


Position 


040 


Fl 
F2 


ΤΟ 


Unit cell ond 





^y 


2 


ኩ 


bith coolant 


En p" 
1.008 1.010 
1.008 1,016 
1.008 1.016 
1.008 1,055 
1.008 1.060 
be 008 1,060 
1.008 VE 
1.008 l.078 
1,008 112: 
1.907 1.190 
1.007 1.251 
1.007 1.320 
1.005 1,381 
1.005 1.1163 
1.005 1.019 
1.005 1.019 
1.005 1.019 
1.00 1.022 
1.968 1.03% 
0.399 1.052 
0,995 1. 100 
7.994 1.128 
0.992 1,170 
0.990 1.800 
0.988 1.270 
1.208 0.388 
1.008 0.997 
1.907 0.982 
1.006 0.981 
1.001, 0,978 
1.003 0.975 
1.003 0.972 
1.000 0.952 
0.3997 0,929 
9.99 0.907 


fu 
1.009 


1.009 
1.000 
1.008 
1.008 
1.008 


1.007 
1.007 
1.006 
1.006 
1,006 


1.005 
1.005 
1.00 


1.007 
1 . 006 
1.001, 
1,004 
1.003 


0.999 
0.995 
0.992 
0.990 
0.959 


0,959 


1,005 
1.90 
1.004 
1.004 
1.003 


1.003 
1,002 
1,000 
0.295 
0.992 


and harmonic eorrection factors 


without coclent 


m” 


k... A 


1.010 


1.042 
1,052 
1.061 
1.061 
1.061 


1.068 
1.072 
1.116 
lue 


1.240 


1.307 
1.368 
1.448 


1.017 
1.020 
1.021 
1.021 
1.032 


1,057 
1,090 
1.107 
1.157 
te TOG 


Le @ipl 


0.995 
0.96), 
0.980 
0.900 
0.576 


23,975 
0,970 
0.952 
0.926 
0.905 
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Table 5 5 Continued) 























Position Te sith coolant Without coolant 
_ ሦ fü p 
HLL 1.001 0,989 0.388 0.988 0.857 
k123 1,001 0,906 0.070 0.985 0,5869 








Ge Veseribtion of Funes in Unit Cell 


in investigating tie flux in the unit cell .uns were 
made along the P, Y and radials emanating from the center 
of the fuel assembly as shown in figure 3. «uns 1 through 
16 were made with no water in the coolant annulus and will 
hereafter be called "dry" runs. „uns 17 throusn 32 were made 
with water in the coolant annulus and will hereafter be 
called "wet" runs. The medium sized foils were used on all 
the dry runs, whereas on the wet runs the foil size was 
varied to study the effect of this parameter on the induced 
activities. The foil spacing was varied on the dry runs but 
was held constant on the wet runs, 

“Spacing 1” is defined as that spacing along a radial 
wien all the foil positions in tne srachite wore Pilied for a 
run. "Spacing 2" corresponded to a foil veing placed in 
every other foil position, and'asaoacinz 3" corresponded to a 


foil being placed in avery third foil position along a given 
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racial. “he above soscing refers only to those folla placed 
in the graphite block. oils were placed in positions in 
the fuel assembly two at a time while the Polls in the 
grachite were being irradiated, and this was called “normal 
spacing” for foils in the fuel assembly. Thus, the data for 
rans 5, 11 and 16 were actually taken during runs 1, 2, 3, 6, 
T, 0, 10, 12, 13 and 1h. The fuel asserbly foil readings 
were grouped inte individual runs almnly for ease of 
reference. At least one run along each racial was mede with 
all or almost all of the foil positions on that radial filled 
both in the fuel assembly and in the graphite. These were 
runs li, 15, 19, 25 and 30, and the foil spacing in the fuel 
asseably on these runs was designated as "close-packed". On 
these runs medium sized indium foils were used. 

Foils ware normally placed along the + and . radials in 
the horizontal position, and along the h radial in the 
vertical position as indicated in Figure 3. On runs 6, & and 
iO medium foils were placed with spacing 1 along the 4 radial 
in a horizontal, a vertical and an leshaped position respec- 
tively. The L-shaped position waa obtained by vending the 
foil into a 90° angle and Inserting it into the olock so that 
1% pointed outward along the radial. 

un runs 27 and 32 along the 2 and kh radials respectively 
the indium foils were wrapved in 0,010-in. cadmium shnest and 


irradiated. 0n these runs only one cadmium wrapped foil was 
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placed in the block at a time in order to aveid too large a 
denression in the thermal neutron flux due to the presence 
of tho cadmium. 

On all but the initial runs the counting timas used 
were either two or three minutes, It was found that there 
was excessive scatter in the experimental dsta using the two 
minute counts, and therefore three minute counts were adopted 
for all the later runs. With tia three minute counts tne 
maximum relative standard deviation in the counting rate vas 
5 per cent with the averate being 3 to ip per cent.  L.xolusive 
of those runs made with close-packed spacing, the foil 
loading for each irradiation averaxed four foils along the 4 
radial and eight foile along the Y and hk radials. Suns were 
made along the C and h radials simultaneously. For any one 
irradiation all the Pcils were counted through once, and then 
a second count was teren., The average of the two saturation 
activities thus obtained was used as a measure of the flux. 
if the foil activities weve high enouch, a third and even 
fourth count vas made and the average of ali saturation 
activities was used, 411 runs meade in the unit cell ave 


listed in Table 6, 
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Table 6. List of runs in unit cell 




















sun Goolant .adial ΡΟΣ} foil "old 
no, size pacing? criontation 
l Hono P medium 1 horízontal 
2 Hone P modium 2 horizontal 
3 Hone P medium 3 horizontal 
iL None P medium 1 ος horizontal 
5 None P dium normal radial 
6 None Q medium 1 horizontal 
? None ህር medium 3 horizontal 
8 Hone 1 medium 1 vertical 
10 None 2 mediun 1 Leshaped 
11 Hone 8) medium normal radial 
Le None i rio dá ura 1 vertical 
13 None t medium 2 vertical 
887 None ig medium 3 vertical 
1S None ፲፻ medium i OF vertical 
16 None R medium Σ radial 
17 «ater P large 2 horizontel 
18 Water Ë medium 2 horizontal 
19 Water P medium l cR horizontal 
20 kater F large/^edium normal radial 
21 water P small normal radial 
22 water Q lar se 3 horizontal 
23 bater à medium 3 horizontal 
eh Water 2 small 3 horizontal 
25 Water d ce dium 1 CP horizontal 
26 Vater Q large/medium normal radial 
27 hater iQ medium horisontal 
28 water R large 2 vertical 
29 Water h me dive: 2 vortical 
30 water Pi medium l CP vertical 
31 water n. large/redium normal radial 
32 water By me dium 2 vertical 








& | 
CP = close packed 
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fable 7.  Hori£ontsi (lux survey at y 7 -10 in., e = 3U in. 











Position x Normalized activity _(c/m) 
(in.) without coolant «ith coolant 
AS -27 205 eli 
BS -21 LOO 512 
CS -15 835 8147 
DS - 9 1088 1062 
ES - 3 1080 1153 
P5 3 121,8 1255 
G5 9 1092 1165 
H5 15 790 869 
I5 21 601 567 
#5 27 zul 268 








D. Horizontal Surveys 


Horizontal vile surveys were mace in the x direetion at 
y = -10 in. and z = 30 in. both with and without water in 
the coolant annulus te determine wiether or not the trans- 
verso flux distribution in the subcritical assembly was 
symmetrical. The normalized activities from these surveys 


are listed in Table 7 and are plotted in migure 29. 
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The raw data for all runs yes reduced to normalized 
saturation activities, Ag, and the corrected activities re- 
ferred to tie center of the uranium sluge Àg’ Tne radial 
distances alons tho P, x and n»n radials were designated p, q 
and r respectively. Various combinations of experimental 
data are plotted in Figures 6 through 29. In fairing curves 
throug the experimental points it was assumed that there 


were no radical changes of curvature of the flux distribution 


within the graphite block. 


A. Effect of Foil Spacing 


“igure 6 indicates tint induced activities for foil 
spncings 1 and 2 along the F redial were approximately the 
same and that they were about 5 per cent lower than the in- 
duced activities oP those folls irradiated at spacing 3. 

This depression increased to approximately 10 per cert for 
those folis located closest to the fuel assembly. along the 
4 radial the change in foil spacing had an effect on the in- 
duced activities as is shown in Figure 7. along the h radial 
there was approximately a S per cent dacrease in induced 


activities of foils irradiated at sgacing 1 compared to those 
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irpadlatod ass spacing 2 as indicated in Figure Ö. 

Un pun | medium foils were close packed on the d radial 
in the fuel assembly and were placed as soacin; 1 in the 
graphite. Figure 9 comosres the activities obtained on this 


c 


in with those obtained en runs l and 5 wecre the foíl 


8pecine 1 was used in tne grapnite and normal spacing was 
used in the fuel asssmuly. There was ap oarontly a 10 to 15 


per cent devression of foli setivity in the fuel esserbiy and 


10 to 20 por cent increase in foil activity in the sraphit 
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The same type of runs was made and cornard on the n radial 
in Fisure 10, There was a 15 to 20 ver cent depression of 
activities in the fuel asserbly with the foile close paced, 
but in the gfaphitbe tie activities were about the same. “he 
above runs were all dry runs. figures 11, 12 and 13 show the 
results of similar runs thit were made with water in the 
coolant annulus. Approximntely a 10 Der cent depression in 
the induced activities was «gasin noted when the foils were 
close packed in the fuel assembly, but there was very little 
chanze in the activitios of thore foils plaoed in the praphite. 


A 


There did not abogar to be a „rent deal of distortion o 


፦ኔ 
ፍተ 
pe 
er 
ሥማ 
Sas 


flux pattern oy placing the file close packed in tho fuel 


assembly. 
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MW. Uifact of Foil Orientation 


kuns 5, 8 and 10 were dry runs made alons the 4 radial 
with medium folles nlased in the graphite in s horizontal 
position, vertical vosition sand an Leshaved position 
respectively. The horizontal and vertical placement gave 
very nearly the sane distribution alors the radial as shown 
in Figura 15. The L-shaped foil orlentition apveared to 
result in activities thet wore depressed aporoximately 10 
per cent from those obtained from the horizontal and vertical 
positions. “ne activities from the Imshased foils also nad a 
larger amount of acatter than those activities obtained from 


foils mounted horizontally or vertically. 


G, úlfect of foil Size 


Large, medium and amall foils were used on runs 22, 23 
and 24 in the graphite on the 2 radial with spacing 3. 
Figure 15 shows that tho small foils resulted in the highest 
specific activity, with the medium foll specific activities 
oeing depressed nporoximately 10 per cent Proe these and the 
large foil activities being depressed 15 to 26 per cent fron 
the mall foil activities. This trend was not observed on 
the P? radisl iere on runs 17 and 18 using lerge and medium 


foils respectively, with spacing ?, almost identical flux 
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distribution curves were obtained along the radial as shown 


in Pigure 16. Figure 17 is a plot of runs 20 and 29 which 


= 


were the sane as runs 17 and 10 except that thoy were taken 
along the H radial. in this instance the large foil 
activities were found to be 12 per cent less than The medium 
foil activities. For run 21 emall foils were placed in tne 
fuel assembly one at a time along the P radial, “ne 
activities obtained in this manner were not appreciably dif- 
ferent from those ootained with medium foils placed two st a 
time in the fuel element assembly, as can oe seon oy oon- 


paring run 21 with run 20 in Figure 16. 


P. Variation of “lux along Different radials 


In an isolated unit cell with cylindrical geormetry the 
lines of conatant flux in the moderator would ve concentric 
circles. In a square unit cell in a reactor the lines of 
constant flux in the moderator in the vicinity of the fuel 
assenoly are closely aporoximated by concentric circles if 
tae overall flux in the reactor were uniform from cell to 
coll. (13, 5» 79) As the unit cell boundary 18 approached 
the lines of constant flux are gradually distorted fron 
circles into souaros. At the cell dboundarythe lines of 
constant flux would be squares. Since the unit cell 


activities, Ags have all been corrested for the cosine 
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distribution in the x direction and the exponential drop in 
the z direction, thase activities correspond to those which 
would be obtained if the flux in the overall assembly were 
uniform. Therefore lines of constant flux (or corrected 
activity) should be very nearly concentric circles near the 
fuel assembly. Sines the lines of constants flux are com- 
pressed slong tho 8 radial as they chance from circular to 
square shape, 15 should be expected that at a given radial 
distance the flux along the P and X radials would be equal 
but less than the flux alone the 4 radial. 

In Figures 15 and 19 the corrected foil activities 
along the F, 4 and K radials are plotted and compared for the 
wet and dry runs respectively. The activities along the dif- 
ferent radials are seen to mateh up very closely for the wet 
runs and fairly well for the dry runs. Activities along the 
Q radial appeared to be slightly higher than on the + and kh 
radials. This was probably due to the reason mentioned above 
and to the fact that spacing 3 was used on the . radial while 
gpacing 2 was used on the + and a radials., The activities 
of the folls in the fuel assembly all fell within relatively 
narrow limits as can be seen in Figures 18 and 19. The 
Single curve faired through these points in the fuel assembly 
corresponds to theo average at a given radial distance of the 
activities measured along the three radials. Figures 20 and 


21 are the same as 1% and 19 exceot that on these runs tho 
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foils in the fuel assembly wors olose pacred, These runs 
simply corroborated the results shown in Figures 18 and 19, 
The average activities in tne fuel assembly are seen to 

agree very closely in general shape, but the magnitudes of 
the close packed foil activities are depressed 10 te 15 per 
cent below the activities of those foils that were irradiated 
only two at a tire. See also Megure 27 which compares the 
averages of the above runs. 
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Ἐν Lfioct of Coolant 


identical runs were made along each radial with and 
without coolant and these are plotted in figures 22, 23 and 
2h. On the wet runs there appeared to be about a 5 per cent 
depression in the flux in the region adjacent to the fuel 
assembly along all three radials. This depression continued 
to the unit cell voundary on the . radial, but there was no 
apparent depression at the unit cell ooundary on the P and kh 
radiaise 

Within the fuel assembly it could in general be said 
that the flux was depressed on the wet runs. Un the + and y 
radials this depression amounted to about 10 per cent whereas 
on the hk radial the depression was very slight. ‘There was a 
characteristic flux pattern evident in the fuel assenbly 


from the data for the individual runs which was more apparent 
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when the average of the readings taren on the three radials 
was plotted as shown in figure 25. There is seen a pattern 
that is fairly consistent with thet predicted theoretically. 
Thet is, for the wet runs there is an increase in thno 
thermal neutron flux across the water annulus, whereas for 
tbe dry runs the flux across the air annulus remains about 
constant. On the P and fi rudials It was noted that on the 
dry runs there wes a rather sharp incroase in flux (about 10 
per cent) in the air hole beyond the process tube, whereas on 
the wet runs there was llttle if any Increase in thermal 
neutron flux aoross this air space. 

Figures 26 and 27 show results of runs which were made 
with and without coolant with foils in the fuel assembly 
close packed and foils in tho graphite at spacing 1. There 
was considerably less scatter in toe experimental data of 
Figures 26 and 27 than there was in Figures 22 and 2h. This 
was probably duo in part to the fact that for the close 
packed runs with spacing 1 in the graphite all of the foils 
were irradiated simultaneously in the same flux field. for 
spacing 2 and 3 there were two or three irradiations required 
to complete a survey along a radial. Since the foils were 
located differently for each irradiation, the flux field was 


snaped differently for each irradiation. 
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F, Cadmium natio 


Cadmium ratios were determined along tne 4 ang i 
radiale with coolant and are plotted in figure 20, along 
both radials the cadmium ratio increased with distance [rom 
the uranium. ‘Tho increase along the „ radial was about 10 
pep cent at the unit cell boundary and along the . radial 
the ratio increased about 17 per cent at the cell boundary 


over what it was in tne fuel aasembly. 


Ge Comparison of Flux in Unit tell With 


Overall “lux in the Assembly 


Table 7 lists the results of the horizontal pile mwrveys 
taken at y = -10 in. and g = 30 in. with and without coolant. 
These surveys vere taken using the large (l in. by 1$ in.) 
aluminum backed indium foils. ¿lots of the horizontal pile 
curveys with and without coolant aprear in Pigure 29. Unit 
cell surveys along the i radial obtained with medium foils at 
spacing 2 are also plotted on Figure 29 to show the relation- 
ship of the unit cell flux distribution to the overall 
assembly flux distribution. The activities per gram obtalned 
with the smaller unit oell foils were apvroximately 1.77 times 
larger than those obtained using the large (1 in. by 13 in.) 
pile survey foils. The unit cell activities were divided oy 


this factor of 1.77, and this reduced unit cell foil activity 
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was plotted on Figure 29. 

Normalized saturation activities in the unit cell along 
the P radial using medium foils at spacing 2 with and without 
coolant were plotted on Figures 4 and 5 respectively to com- 
pare the vertical flux distribution in the unit cell with 
that in the pile. Since the overall vertical pile survey was 
taken at x = -3 in. and the vertical unit cell survey was 
taken at x = -10 in., there were no common points in the two 
surveys by which one set of data coulá be converted to the 
other. It was assumed that the ratio of activities would be 
the seme as it was for the horizontal surveys. Thus the unit 
cell activities plotted in Figures lh and 5 are the actual 
saturation activities reduced by the factor of 1,77 and 
corrected for the cosine distribution in the pile. The 
eorrected vertical flux survey in the unit cell matched the 
vertical pile survey in the region remote from tha fuel 


assembly. 
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VIII. DISCUSSION OF AESULTS 


A. Techniques for Measuring Flux Distribution 


in a Unit Cell 


In making any physical measursment great caro must be 
taken not to influence the quantity measured by the technique 
of measuring. This is especially true in the case of neutron 
flux measurement by the foil activation method. It is im- 
possible to avoid altering the flux field when using the 
activation method, out it is possible to keen these altora- 
tions as small asa possible. Space must be made available for 
placing the indium foils, and the rolls must be large enough 
ao that their induced activities after irradiation in the 
neutron flux are detectavle and meaningful. 

in order to obtain reproducible results it is of primary 
importance to odtain sufficiently geod counting statistics. 
Due to the smaliness of the indium foils, the low neutron 
flux, and the 5h minute half Life of indium, the allowed 
counting time is limited. Use of gold foil with its longer 
half life would help solve this problem, but would require 
Longer irradiation times. If several foils were to be 
counted it was found best to use short counts and count 


through all the foils two or more tires and average these 





rather than take single long counts. although averaging the 
saturation activities obtained from several short counts did 
not appreciably improve the counting statistics it did help a 
great deal in reducing the scatter of the experimental data. 
This was probably due to the fact that averaging tended to 
reduce the variations in countins geometry from one foil to 
another. The countins gewmetry of the foils must of course 
be kept constant as well as the geometry of the foll while it 
ig being irradiated., This was probably tae main reason why 
the L-shaped oils gave such erratic results on run 10. in 
bending the foils into the 50° angle some were undoubtedly 
bent slightly different than others and thus nad different 
irradiation geometry. Furthermore when these tent foils were 
counted it was difficult to flatten thom out under the 
counter and this introduced variations in counting geometry. 
By irradiating the foils flat, either horizontally or 
vertically, the above difficulties were eliminated. 

Foll size, spacing, loading and counting time should be 
determined by the type of results desired ang the time 
available to collect the data. If scuantitative resuits on 
the unit cell are desired, the smallest foils possible should 
be used, and they should se irradiated in the unit cell ons 
at a tine. The minimum size of the foil would depend on the 
flux level. For the flux encounteced in these experiments 


it vas ound that tho minimum usable size of U.003-in. indium 
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foil was ë im, by + in. If only cualitative results aro de- 
sired, for example, if it is desired to determine only the 
pattern of the flux distribution, the foils placed in the 
unit cell may be larger both in size and in number» Even 
piscing the foils close packed in the fuel assembly and 
placing the foils at spacing 1 in the graphite cid not appear 
to greatly distort the flux distribution in the unit cell. 
Althougn the general flux level was depressed by close packing 
the foils in tae fuel assembly, the general shape of the flux 
distribution was not greatly altered as is shown in figure 
25. 

Using the activation method is a time consuming process 
wnen using small foils in a low neutron fiux due to the fact 
that the foils must be irradiated about 6 hours vetween each 
run. It is therefore most advantageous to read as many foils 
as possible per irradiation. With the size oiis used it was 
found that the maximr numbor of feila it was practical to 
read at one time was elyht. By using two or moro counters 
símultancously it would be possible to reduce the time re- 
quired for counting per irradiation and also the number of 


irradiations per run, 


B, The Flux Distribution in a Unit Cell 


The matching of tho flux distribution curves along the 
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de” 


s Qand H radials as was seen in Pírures 18 and 19 would 
seem to verify the validity of the f, and f, correction 
factors which were questioned by Clayton (2, ». 33). At 
least they seem to ne valid within the limits of accuracy of 
this investigation. 

The ganeral effect of water in tho coolant annulus of 
the fuel &ssenbly was to depress the flux in the fuel assembly 
and in the surrounding graphite. With the exception of the 
flux pattern shown in Figure 23 it apoesarod that there was 
very little if any devression of flux in the region near the 
unit cell boundary. Hence, horizontal and vertical flux sur- 
veys made across the assembly with the survey foils located 
Near the unit coll boundáeries would not detect any variation 
in tho flux distribution due to the addition of coolant., The 
survey foil positions in the subcritical assembly used in 
this investigation were located half way between the centers 
and the boundaries of the unit cells so tnat the depression 
of flux due to the coolant was hardly detectable. although 
not very pronounced tho moderating effect of tha water was 
apparent from the experimental data plotted in Figures 23 
through 26. With 1noroved statistics and refined counting 
procedures it should be possible to measure quite ascurately 
the effect of coolant on the flux distrioution in the fuel 
assembly. 


The hish foil activities obtained in the unit cell as 
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comoapad with those activitics obtsined when making hori- 
sental flux surveys across the entire assenbly wore probably 
due mainly to the counting geometry. The geometry factor Por 
the unit cell foils aprreached 50 pex cont since these foils 
were completely covered by the window of the counting tube 
which was 1 1/8 in. in dianeter. However the large foils 
used in making the overall pile surveys extended out beyond 
the counter window. This would cause the activity per gran 
for the small foils to be higher than for the largo roils. 
The unit cell foil activities which ere plotted on 
Figuro 29 indicate that the flux distribution across the 
asambiy has large deviations from the cosine distrioution 
due to the depressions ln the vicinity of tne fuel assemolies. 
This points us again the importance of foil placement when 
making flux surveys. Three different horizontal flux dis- 
tributions would be obtained across the agsambly depending on 
whether the foils were placed in the emoty holes, in the sur- 
vey slots or in thse holes loaded with fuel ele=ants, These 
are shown as curves à, B and C respectively in Figure 30. 
The actual flux distribution would be shaved lize curve ሀ 
with depressions at the fuel olements and peaks at the cell 
boundaries. The depressions in the region of the fuel 
assemblies are probably not as pronounced as the experimental 
data indicates because the presence of the indium was partly 


responsible for the lowered flux. 
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Unit cell foil activities along the P radial were plotted 
on Piígures l| and 5 simply to show that the vertical flux dis- 
trioution also deviztes considorably [rom the theo:stical ex- 
ponential drop, It is apparent that when irradiating foils 
for the purpose ef determining the buckling, it is important 
that the foils be irradiated at the same relative position in 
cach unit cell. The validity of the position correction 
factor, -— would again appear tc bo verified from tne fact 
that when it was applied to tne unit cell survey data, the 
corrected data matched the pile survey data very closely in 
the region away fron tho fuol assombiy. 

The shapes of the cadmium ratio curves in Pimire 28 show 
that the fast neutron flux in the region near the fuel 
assembly is higher than in the region near the unit cell 
boundary. The cadmium ratios obtained along the < radial 
indicate that the fast nevtron flux is syrrretrical with 
respect to the center of the fuel slug. Since a low cadmium 
ratio indicates a relatively higher fast neutron flux, it is 
apparent that the fast flux becomes maximum in the fuel 
assembly, dro»s off os tho col' boundary is aporoached, and 
increases again as the slug in the adjacent unit cell is 
approached. The higher cadmium ratio at the cell boundary 
along the i radial as compared to that along the radial wes 
due to the factthat the cell center-to-boundary distance was 


greater along the R radial, sand hence there were fewer fast 
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neutrors remalninz nt the cell boundary in the r. direction 


than there were in the „ direetion. 


C. Comparison of „xperimenteal .esults with Theory 


Tue theoretical flux distribution in the unit cell of a 
simple two- region fuel-moderator system is shown as curve A 
in Figures 31 and 32. Curves E and C are the theoretical 
flux distributions in a multiregion system with and without 
coolant respectively. Since the exverimental results indi- 
cated there was very little depression of flux in the region 
near the unit cell boundary due to the water coolant, the 
theoretical curvos were normalized to the coll bounda:y flux. 

::989ፌ 1508 of flux surveys along the PF and hk radials for 
both wet and dry runs are plotted on Miyures 31 ard 32, where 
the experimental data was normalized to the average ceil 
boundary flux. Ywoeregion theory is seen to agree quito well 
with the experimental data. However a curve throusn the 
experimental points wouid have less slope than that predicted 
by kurray. Multiregion theory both with and without water 
coolant appears to give hich values for tne flux in the 
rrepnite bub the shane of the curve appears to agree closely 
with the experimental data. In multirepion theory the varia- 
tion of flux across the cladding, coolant ané process tudes 


wes assumed to be linear, and tuo flux across air gaps was 
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Figure 32. Comparison of experimental data with theory 





10%. 


assumed to be constant. These assunntions eapyear to ba valid 
except for the water coolant. The theoretical decrease in 
flux across the water annulus was about twice as large as 
that observed. This discrepancy was prooaoily due to tne fact 
that the moderating effect of the water was ignored in 
Calculating the theoretical curves. Although the theory of 
humsey and Volkoff provides for tho moderating effect of the 
water in determining the thermal utilisation of the lattice, 
it does not provide for the inclusion of this effect in 


calculating the flux distribution. 





JX. CONCLUSIONS 


The activation method is suite adastable to the measure- 
ment of neutron flux in tne unit cell. Sy reducing the 
statistical deviation of the foil activities it should oe 
possible to determine accurately the flux diatribution. For 
more exact calculations it would also be necessary to correct 
foil activities for the fast neutron flux which is present 
and measured along with the thermal flux. 

The flux distribution predicted by two-region theory is 
in very close agreement with the observed districution. 
Murray's assumption that any poisons thet are tolerable go not 
appreciably disturb the basic fuel-moderator flux distribution 
appears to be reasonably valid, although these poisens de 
flatten out the flux distrivution in the graphite. Piti- 
region theory appears to predict a flux level in the 
moderator which is generally higher than was observed. It 
also gave a larger flux depression across the water coolant, 


because the moderating effect of’ the water was ignored. 
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further investi;ation of the flux distribution in tne 
it coll of the suceritical Asgerbly could te carried out by 

varying certsin osher parameters, such as coolant, lattice 
size and slug size. Une could also examine otner unit cells 
either adjacent to the ons examined herein or in unother 
region of the assembly. ùin any further work tno flux surveys 
might be continued into the fuel elerent itself. 

Tne aprarent rise in the flux in tue air zao beyond the 
process tuse along the í; and lt radials presenta an interesting 
phenomenon which could be further investirated. The flux 


pattern in this region is apparently dexsendant upon the flux 


k 


a 


pattern in the coolant amnulus, however the statistical de- 

viation of the experimental data of this investigation sro- 

hibits making any definite conelusion alorsz tnese lines. 

ánother suojcot fer investigation could be toe theoretical 

evelopment of zhe flux distribution in the unit gell. As 

was previously cointed cut, the theory of iunsey and Volxoff 

was primarily aimed at a more refined prediction of the 

taermal utilisation rather than an exact solution of the 

point-to-point flux distribution. By extending diffusion 


theory to the multiresion system it would be possible to 
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obtain the theoretical flux distribution in each region, 
Such treatment would be particularly adaotable to precicting 
the flux distribution scross a moderating region, such as 


» > 


the water filled coolant annulus, 
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Table 06, Dimensions and matoriel constante for the unit cell 

















Dinmensiona 

Ty 9 uranium rod radius 1.270 cm 

tal s bhickness of alurinum cladding 0.202 ατα 

Uu > effective thickness of water annulus 0.273 cm 

to ; 6ffoctive thickness of vrocess tube 0,102 cm 

tagps effective thickness of air annulus 0.1455 cm 

ry ፣ equivalent inner radius of graphite Cod ar 

Po >, Oduivnlent outer radius of graphite 11.91, om 
Volume per siug 

Va» uranium 103.0 om? 

Maz » Slur eAn and cap 28.) cr 

ú ‚water 5557 cm 

Vo » procesa tubo 22,0 cn? 

Es . £rap»hito 9250 cm? 
Aosorotion cross sections 

Sa,» aluminum 0.01323 cm”? 

S a uranium 0. 32h om”? 

Ze , graphite 0.00036 cm”? 

sie water 0.017 di 
Inverse diffusion lLenrths 

Ky» uranium 0.675 cm”! 

K . » craphite 0.01992 cm”? 

X u s water OZ enr i 


X a3? aluminum 0.0195 omil 
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